Thermodynamic Properties of the SU(2)f Chiral 

Quark-Loop Soliton 

M. Schleif and R. Wiinsch* 

Institut fiir Kern- und Hadronenphysik, Forschungszentrum Rossendorf e.V., 
Postfach 51 01 19, D-01314 Dresden, Germany 



Abstract 

We consider a chiral one-loop hedgehog soliton of the bosonized SU(2)f 
Nambu & Jona-Lasinio model which is embedded in a hot medium of 
constituent quarks. Energy and radius of the sohton are determined in 
self-consistent mean-field approximation. Quasi-classical corrections to 
the soliton energy are derived by means of the pushing and cranking 
approaches. The corresponding inertial parameters are evaluated. It is 
shown that the inertial mass is equivalent to the total internal energy of 
the soliton. Corrected nucleon and A isobar masses are calculated in de- 
pendence on temperature and density of the medium. As a result of the 
self-consistently determined internal structure of the soliton the scaling 
between constituent quark mass, soliton mass and radius is noticeably dis- 
turbed. 

PACS: 12.39.Fe, 12.38.Lg, 12.40.Yx, 24.85.+p 
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1 Introduction 



Chiral soliton models have proven to be a fruitful approach to the description of 
nucleon structure. Starting from the Nambu & Jona-Lasinio (NJL) lagrangian 
[|l| and applying a well defined scheme of approximations one was able to obtain 
stationary and localized field configurations denoted as non-topological chiral one- 
loop NJL solitons. They can be used to model nucleons, A isobars and strange 
baryons on the basis of interacting quarks (for review see [@, |^). 

The NJL lagrangian incorporates chiral symmetry and its spontaneous break- 
down IQ. It has been used to study the restoration of chiral symmetry in a hot 
and dense nuclear medium modeled by a gas of constituent quarks (for review see 
0). The decrease of the constituent quark mass at higher temperature and/or 
density of the medium describes the phase transition from the chiral condensate 
to the chirally symmetric phase. The calculated effects are in satisfactory agree- 
ment with the predictions of lattice calculations and of the chiral perturbation 
theory as well. 

It is an attractive idea to combine both features of the NJL lagrangian and 
to study the behavior of a soliton embedded in a hot gas of constituent quarks 
with a dynamically generated mass. Such a model incorporates the restoration 
of chiral symmetry and the possible dissolution of the soliton, which simulates 
the deconfinement transition of hadronic matter. In contrast to many other ap- 
proaches studying medium modifications the non-topological soliton model equips 
the baryon with an internal structure which may be modified by the medium. 

Using this approach as a model for baryons in hot hadronic matter one should 
be aware of its approximative character which is even not free of inconsistencies. 
Below the critical values of temperature and density, the quark gas is not the 
ground state of strongly interacting matter, neither in nature nor within the 
model. If the soliton is stable the medium itself consists of solitons. This goes 
beyond the mean-field approach. The effect we can study within a mean-field 
picture is the scale change connected with the reduction of the constituent quark 
mass at increasing values of temperature and density and its effect on the self- 
consistent mean-field. Such an approach rests on the assumption that the dom- 
inating effect of the medium consists in the reduction of the constituent quark 
mass while its local variation is of minor importance. The free motion of the 
quarks representing the medium as a quark gas is an obvious shortcoming of 
the approach and may overestimate the influence of the medium on the soliton. 
There are attempts ||^, to replace the quark degrees of freedom in a part of the 
effective action by nucleonic ones without introducing new parameters. The re- 
sults are not very encouraging since chiral symmetry is restored already at normal 
nuclear density in this approach For a more detailed discussion see Ref. [§. 

The soliton which we investigate is in most respects identical with the soliton 
described in Refs. [Q, |^, ^. The differences concern the particular treatment of 
the valence quark level and the use of the chemical potential for adjusting the 
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baryon number of the soliton. 

Due to mean-field approximation and hedgefiog ansatz the soliton has de- 
fects already known from the soliton in vacuum: it violates translational and 
(iso-) rotational invariance. Therefore it is affected by center-of-mass motion and 
represents a mixture of nucleon and A isobar instead of a particle with definite 
spin and isospin. The violated translational and rotational symmetries can ap- 
proximately be restored. The quasi-classical pushing and cranking approaches 
1^ constitute a feasible way to exclude spurious contributions to the energy and 



to equip the soliton with the correct values of spin and isospin. The size of push- 
ing and cranking corrections is controlled by inertial parameters. While we relate 
the inertial soliton mass to its total mean-field energy the (iso-) rotational moment 
of inertia is calculated numerically. The relation between inertial mass and inter- 
nal energy, which is derived in this paper, is an extension of the corresponding 



relation for a soliton in vacuum IITT 



In Sect.|^, we shortly outline the basic ideas defining the NJL soliton in a 
medium of constituent quarks at finite temperature and density and review the 
main formulae. We determine that region of density and temperature where a 
stable soliton exists. The baryon number of the soliton and its spatial distribution 
is considered in Sect. 3. Here we critically discuss the method to fix the baryon 
number to one, which was applied in Ref. The numerically determined soliton 
energies and radii are given and discussed in Sect. 4. In Sect. 5, we determine 
quasi-classical corrections to the soliton energy. We consider the soliton in a 
boosted and rotating frame and calculate the corresponding inertial parameters 
and energies. The corrected nucleon energies are given and discussed in Sect. 6. 
Conclusion are drawn in Sect. 7. An appendix completes the calculations in 
Sects. 2, 3 and 5. 



2 NJL soliton in a heat bath 

We consider an ensemble of up and down quarks with Nc = 3 colors and an av- 
erage current mass m = [m^ + at temperature T and chemical potential 
/i = /iu = /id- The latter will be related to temperature T and density po of the 
medium embedding the soliton. The quarks interact via a four-quark contact in- 
teraction, which consists of a chirally symmetric combination of a scalar-isoscalar 
and a pseudoscalar-isovector term, with the coupling strength G/2 introduced by 
Nambu & Jona-Lasinio The soliton is defined by an effective action whose 
derivation from the SU(2)f NJL Lagrangian incorporates the following steps (for 
a review see Refs. 0, 0]): 

1. Introduction of auxiliary meson fields a and tt by means of a Hubbard- 
Stratonovich transformation O] in the generating functional using the 



imaginary-time formalism. 



2 



2. Derivation of an effective meson action ^cfr[o", tt] by applying the stationary 
phase approximation on the meson fields (no meson loops, a and tt as 
classical mean fields). The effective action obtained in this way consists of 
a purely mesonic part and of a fermionic part A'^. The latter describes 
the contributions of the various quark levels to the effective action (quark 
determinant). 

3. Restriction of the meson fields to static and spherically symmetric hedge- 
hog configurations (cr(r,r) = cr(r), 7r(r,r) = 7r(r)f). In our numerical 
calculations, the meson fields will additionally be restricted to the chiral 
circle (o"^(r) + 7r^(r) = cig = const). Otherwise, a stable soliton does not 
exist 

4. Splitting the quark part of the effective action into a contribution A'^'^'^'^ 
{sea contribution) which results from a completely occupied Dirac sea and 
a residual contribution A^''^^'^(T, fi) {medium contribution) which describes 
the occupation of the quark levels according to temperature and chem- 
ical potential (quarks in levels with positive energy and holes at nega- 
tive energy). The sea contribution diverges and is regulated by means 



of Schwinger's proper-time regularization scheme [16|. The corresponding 
cut-off A is not considered as a free parameter but is related to the exper- 
imental values of the pion mass and of the weak pion-decay constant in 
vacuum B, \i7\, ^ 



5. Interaction strength G and cut-off parameter A are determined in the vac- 
uum and assumed to be independent of T and /i. This assumption ensures 
the exact scaling between pion decay constant and constituent quark mass. 

6. The soliton itself is defined as a localized deviation of the fields from their 
asymptotic values ctq and ttq = which describe the homogeneous medium. 
Solitonic expectation values are defined by the difference between the values 
obtained for solitonic and homogeneous field configurations. 

The relative effective action of the soliton is obtained by subtracting the effective 
action ^efrfcojO] of the homogeneous configuration from the effective action of 
the solitonic field 

Acs[cr, tt; o-q] = Aff[o", tt] - Acs[cro, 0] = ^"[cr, tt; o-q] + A'^W, n; (Tq] • (1) 
It consists of a purely mesonic part 

A"'[a, ^o] = / dV [a'{r) + n\r) - a^] +^^J dV K - a{r)] (2) 

and of the quark determinant which can be written 

^^[a,7r;ao](T,/i) = -iV,Trln^^ = ^^'^^^[a, tt; ao] +^^''"^"[a, tt; ao](T, /x) (3) 
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with 

A^'^^^[a,n;ao] = "^^c lim TTr In ^ (4) 



and 



^'^•'"^'^[fT,7r;ao](T,/i) = - ^'''^^^ (5) 

- -iV.Trln^ + iiV.limTTrln.^^°) 



with the trace Tr defined in appendix A. While the medium contribution (|^) 
is finite and vanishes in the hmit (T, /i) — >■ the sea contribution (^) diverges 
and does not exphcitly depend on the thermodynamical variables. The latter is 
regularized by replacing the operator trace Tr ( [A .21 ) by a regularized trace TrA 
(|A.3|) . The single-particle operators 

= dr + h-fi, (6) 
Doi^i) = dr + ho-fi (7) 

consist of the derivative dr with respect to the euclidean time coordinate r, the 
quark hamiltonians 

h = h{cx,TT) = a-p + /3 [(T(r) + i75T-f 7r(r)] , (8) 
ho = /i(cxo,0) = a-p + (3ao, (9) 

and the chemical potential fi. The Dirac matrices are denoted by /5 = 7°, 7 = 
(7"'^, 7^, 7'^), 75 = i7°7^7^7^, cx = Pj, and r is the vector of Pauh matrices. 
Spatial coordinates are denoted by r and have the components r\ the absolute 
value r = |r| and unit vector r = r/r. 

The crucial quantity for the description of a grand canonical ensemble of 
quarks is the thermodynamical (grand canonical) potential given by 

fi(T, jj) = TAff = n'^ + n'^iT, /i) (10) 

with 

On the analogy of the effective action we split the quark part of the canonical 
potential into and a medium contribution 

Q'^iT, fi) = -N.TTta In = VlY"" + n^'"^^<^(T, /i) (12) 

where TrA means regularization of only the sea contribution. For time-independent 
meson fields the determinants of the inverse propagators (^ ^ are real and the 
regularized sea contribution can be written 

f^r-^-^HmTTrAln^M^). (13) 
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In the proper-time scheme, we get by means of Eq. ( |A.3| ) 



(-)q,sea 
"A 



Nc f ds r auj y-^ r 
~9' . ~7 . 97^1 



1/A2 -oo 



-^E A) \ea\ - i?E(e°, A) 



(14) 



where (5°) are the eigenvalues of the quark hamiltonians h [ho) defined in 
Eqs. (H, H), and i?E is the regularization function 



i?E(£,A) 



1 £^ 

'2' A2 



(15) 



with the incomplete Gammafunction r(x, a). Notice that the degeneration with 
respect to the color degree of freedom is explicitly taken into account by the 
factor Nc and included neither in the trace Tr nor in the sum over a. 

The medium contribution to the quark part of the canonical potential (|12D is 
finite and will not be regularized. One gets by means of Eqs. (|A.lHA.5l) 



-AT.TTrln , , 



+ Nc lim TTrln 



/^o(0) 



1 _|_ p-sign(eQ) {ea-^i)/T 

-Nc^i B"""" - NcT V In -— . 



(16) 



The medium contribution depends on the thermal occupation probability of the 
various quark levels which are controlled by temperature and chemical potential. 
The quantity 

^sea^_^sign(..) (17) 



describes the baryon number of the Dirac sea for the solitonic field. Usually 
the number of quark levels with positive and negative energy are equal and B'^^'^ 
vanishes. It differs from zero only if the meson field is strong enough to pull down 
one or more quark levels from the positive continuum into the negative energy 
region. This happens at rather large interaction strength G corresponding to 
vacuum constituent quark masses M > 700 MeV, and we shall not consider this 
case here. 

Customarily one treats the contribution of the valence level (Q;=val) to the 
medium part ( |T6| ) separately, ascribes occupation number one to this level {n^^^^ = 
1) and leaves it empty in the homogeneous medium (n^o ^ = 0) 0. This is the 
simplest way to realize a soliton with baryon number one in a cold medium. 
However, the hole in the homogeneous configuration has serious consequences for 
the size of the iso-rotational moment of inertia which will be studied in Sect.lF^. 
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The rule to regularize only the sea contribution to the quark determinant 
should be considered as an ingredient of the model. It does not reproduce the 
correct limit T — > cxd but dealing with a low-energy model we need not consider 
this case. In our case, the regularization procedure would have a negligible effect 
on the medium contribution since the cut-off is larger than chemical potential 
and temperature (A > fi + T). Moreover it simplifies the model considerably 
since it decouples the regularization procedure from temperature and density 
dependence. 

The classical meson fields a and vr minimize the grand canonical potential 

m 

r ( ^ = and = 18 

leading to the equations of motion 

a{r) = m-G{{q{r)q{r))), (19) 
7r{r) = -G{{q{r)h,r-rq{r))) . (20) 

In general, the equations of motion can only numerically be solved since the 
thermal expectation values ((...)) on the right sides depend functionally on the 
fields on the left sides. Expectation values of currents such as in Eqs. (|T9|, ^) 
will be evaluated in Sect.^. A particular solution of the equations of motion is 
given by homogeneous fields a(r) = (Tq and 7r(r) =0 where uo has to fulfill the 
gap equation which follows from Eq. (|19]). A constant sigma field acts as a mass 
on the quarks and ao{T,^) is identified with the constituent quark mass M* . Its 
value M at T = /i = is the only free parameter of the model, which can vary 
within reasonable limits (see e. g. 0). It determines the strength G of the quark- 
quark interaction in the initial NJL lagrangian. Keeping G fixed the constituent 
mass M* for finite values of temperature and density is uniquely determined by 
the gap equation. We chose M = 420 MeV in the numerical calculations. This 
value reproduces the experimental A-nucleon splitting. 

A solution of the equation of motion is called a self-consistent field configura- 
tion since one considers not only the explicit dependence of Q (|T0|) on the meson 
fields via fi™ but also the dependence via energy spectrum {sa} of the quarks 
which enters the parts Q^^'^^ (^) and fi*^'™°'^ (0)- Restricting the meson fields 
to the chiral circle a and vr fields are not independent of each other and equa- 
tions (p!9| , pop can be replaced by a single one e. g. for the chiral angle d{r) (see 
e.g. ||19[). We consider hedgehog fields with winding number one characterized 
by the boundary conditions 9(r = 0) = —n and 9{r—^ oo) = 0. 

The lack of confinement in the NJL model forces us to exclude the valence 
level from the thermal equilibrium and to keep its occupation probability fixed to 
one independently of temperature and chemical potential as proposed in Ref. 0. 
The valence quarks play a crucial role for the existence of self-consistent solitonic 
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field configurations. Only the valence quarks yield a spatially restricted negative 
contribution to the expectation value on the right side of the equation of motion 
(pJl) leading to a well in the a field. The soliton is stable if the well is deep 
enough to bind the valence quarks. If one occupies the valence level according 
to the thermal occupation probability, which is smaller than one, the resulting 
well binds the quarks weaker, and - starting from a critical temperature - a 
homogeneous field with free quarks is the only self-consistent solution of the 
equations of motion. This happens already at temperatures around 100 MeV 
far away from the expected transition point to the quark plasma. Keeping the 
occupation number of the valence level fixed the plasma transition takes place at 
reasonable temperatures around 180 MeV . This transition does not coincide with 
the restoration of chiral symmetry indicated by the reduction of the constituent 
quark mass M* to the value of the current mass m. The constituent mass is only 
reduced to half of its vacuum value when the soliton dissolves. 

Fig. 1 outlines that region in the T—pQ plane where we have obtained stable, 
self-consistent solitonic field configurations. The medium density po is related 
to T and fi via Eq. (|45|). The region with temperatures T < 75 MeV (below the 
broken line in Fig. 1) has to be considered with some caution since we performed 
our numerical calculations within a discrete basis pO[ by introducing a box with 
radius D. Below 75 MeV, the meson fields start to oscillate during the iteration 
and the final results are very sensitive to the box radius. The finite box radius 
produces an artificial spacing and shift of the quark levels which are proportional 
to 1/D. Shift and spacing are important in a transition region around the Fermi 
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Fig. 1: Region in the 
T — Po plane where solitonic 
field configurations have been 
found for M = 420 MeV. The 
density po is given in units 
of the normal nuclear density 
/Onm = 0.16fm~^. In the region 
below the broken line the self- 
consistent meson fields exhibit 
pronounced oscillations out- 
side the soliton in the course 
of the iterative solution of the 
equation of motion. 



energy where the occupation probability varies rapidly. The width of the tran- 
sition region is proportional to the temperature. In the course of the iteration 
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the levels in the sensitive transition region around the Fermi level change rapidly 
their contribution to the mean field with a significant effect on its shape. The 
calculation is stable if a larger number of levels lies within the transition region, 
i. e. if the level spacing is sufficiently smaller than the transition region. At low 
temperatures the spacing has to be rather small and the basis for a reliable cal- 
culation must be large. In this way the capacity of the computer determines a 
lower temperature limit for a reliable calculation. We used a box with radius 
D = 18/M* which restricts ourselves to temperatures above the broken line in 
Fig. 1. In contrast to finite medium density a calculation at vanishing density is 
not affected by the level spacing. In this case, the Fermi energy lies in the middle 
of the energy gap between ±M* and there are no quark levels in the sensitive 
region. 

At temperatures and densities above the solid line in Fig. 1 a solitonic solution 
of the equations of motion (|I^, ^) has not been found. Here the self-consistent 
meson field is too shallow to bind quarks. 

Knowing the grand canonical potential Q the free energy F of the soliton 
can be obtained by means of a Legendre transformations replacing the inde- 
pendent variable fi by the baryon number B = —dQ/{Ncdfi). The internal 
energy is obtained by an additional Legendre transformation from the depen- 
dence on temperature to entropy S = —dfl/dT. Analogously to the effective 
action we split internal and free energy into mesonic, regularized quark-sea and 
quark-medium contributions and subtract the corresponding energies of the ho- 
mogeneous medium 



E 
F 



E™ + EY^"^ + E"^'™""^ 



(21) 
(22) 



Since mesonic and sea contributions to the grand canonical potential are inde- 
pendent of T and n we have 



and 



7-iq,sea T-iq,sea Qq,sea 

-^A — -^A — "A 



The medium contributions are given by 

^q,med _ 



^q,med^y^^^ 



(23) 
(24) 

(25) 



and 



med 



1 

on 



(26) 
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a 

-N^Yl [^ln(l-sign(£°)n^o(r,/i)) +^n^o(T,/i) 
where we have introduced the modified occupation number 

~ ,rj. ^ ^ r.(_ \ _ Sign(£a) 

which describes the thermodynamical probabihty to find an occupied level at 
positive energy 8^ and a hole at negative energy, respectively. The latter is 
supplied with a minus sign. For the completely occupied Dirac sea without any 
additional quarks above it we have ^^^(0, 0) = Va. 



3 Baryon number, density and chemical poten- 
tial 

Now let us investigate the baryon number B of the self-consistently determined 
solitonic field and their spatial distribution p{r). For that aim we consider ther- 
mal expectation values ((O)) of one-body quark operators 



= jd^rq\r)Oq{r) {2i 



where (9 is a time-independent operator acting in the Dirac and/or fiavor (isospin) 
space. The baryon number is obtained with = 1/Nc. 

To calculate thermal expectation values of one-body quark operators (^) we 
define a generating function 

fi?^)(T,/x;K) = -iV,TTr(A) In = + fi^'-^^(T, /i; «:) (29) 

given by the canonical quark potential (0) with the inverse propagators D(o) (yu) 
replaced by 

L'(o) (/i; k) = D^o) (fi) - kO = dr + /i(o) - /i - . (30) 

Restricting the meson fields to their classical values the mesonic part of the grand 
canonical potential does not infiuence expectation values. We shall use both the 
unregularized version 
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and the regularized version Q^'^^i'^) of the sea contribution to ( pQ]) with Tr re- 
placed by TrA. The medium contribution to the extended canonical potential 



(29) is given by 



-iV^TTr In + lim TTr In 



(32) 



Expectation values of an operator (p8|) can be expressed by 



K=0 



mod 



(33) 



with the unregularized sea contribution 



{or'' = - 



K=0 



—Nc lim TTr 



DiO)-'-DoiO)-')o] (34) 



- V ^ C'a-sign(e°) 



and the medium contribution 

,^ed _ dn^'^'''{T,^I■,K) 



my 



dn 



(35) 



K=0 



= iV,5:[n,jr,/i)0«-n,o(T,/i)0^ 

with the modified occupation numbers n (o) (T, /i) (p^) and the matrix elements 

= («(o) |o|«(o)^ = |dV$i°)t(^)o$(o)(^) (36) 

of the operator O with the normalized eigenfunctions ^^^^{r) of the hamiltonian 



h (ho). Sea contributions such as expression (|3^) are defined as expectation values 
at zero temperature and we use the single brackets instead of the double ones 
which stand for a thermal expectation value. In fact, the sea contribution is not 
completely independent of T and but depends on them via the self-consistent 
mean fields a and vr. Using the regularized version of the sea contribution (pl|) 
we get the regularized expectation value 



{O) 



A 



(37) 



In the proper-time scheme, the regularization function is given by 

sign(e) ^ /l £2 \ r / / A N 
i?m = ^7^r -,- =erfe(./A) 



(3^ 
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with the complementary error-function erfc(x) = ^ J^dte Inserting O = 

1/Nc one gets the solitonic baryon number 

^ = {{W f ^(^))) = + ^ - (39) 



with the unregularized sea contribution introduced in Eq. The same expres- 
sion is obtained if one starts from the grand canonical potential (p!OD and uses 
the thermodynamical relation B = —dQ / {N^ dfi) keeping in mind that the meson 
fields have to minimize the potential (pISl). 

To investigate the properties of a soliton which is embedded in a medium with 
given density po we have to establish a relation between T, po and fi. This will 



be done below (45). Knowing T and /i one can determine the solitonic field by 
means of the equations of motion (|IU], |2y). Its baryon number (^) varies with T 
and /i and is different from one in general. The usual method to get a state with 
definite baryon number by minimizing the free energy can not be applied since it 
changes the chemical potential which has already uniquely been determined by 
the medium density po- In Ref. 0, a chemical potential /ig for the solitonic field 
configuration was introduced, which differs from the chemical potential /x of the 
homogeneous field, in order to fix the solitonic baryon number exactly to one. 
However, such a soliton is spatially unlimited since a finite fraction of the baryon 
number is uniformly spread over the whole space. To elucidate this statement 
we consider the baryon density which is defined as the expectation value of the 
current 

0{r) = q^{r)Oq{r). (40) 

with O = l/Nc- The expectation value of currents (1^) with a time- indepen- 
dent operator O can be treated in a way similar to the expectation value of the 
operator (pHj). One defines a generating functional Q'^[k\{T,p) by formally the 
same expression ( ^9]) but with a space-dependent function n{r) instead of the 
parameter k. The corresponding expectation values are obtained by Eqs. (pl|-p8|) 
with the derivative d/dn replaced by the functional derivative 6/6K,{r) and the 
matrix elements 

0^^\r) = ^^^^\r)0^^^\r) (41) 

instead of the matrix elements (|5B|). The expectation values in the equations 
of motion (0, [2^) are of the same type and can be obtained with C = 7° and 
O = i7°75r-f, respectively. Applied to the baryon density we get 

l/T 

p{r) = -T ydr(rr|tr[D(p)-^-Do(Ai)"^]|^'^) =p''''(r')+p"'''^(r*) (42) 



with 

^sea^^^ = -V[sign(eJ<|.t(r)<|.„(r)-sign(5°)<|.°t(r)$0(r)], (43) 
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^med(^) = E [neAT.^^) ^^r) $.(r) - n,o (T,/i) $ot(r) . (44) 



Integrating over the whole space we recover the total baryon number (pQf ). 

First let us consider the homogeneous medium characterized by the hamilto- 
nian with a constant a field o"o = M* and vanishing vr field. The corresponding 
eigenfunctions are plane waves characterized by the momentum vector k and 
normalized to one particle in the volume V. The sea contribution (|43|) vanishes, 
and the sum in the medium contribution (^4|) has to be replaced by an integral 



4V / taking into account both signs of the energies ie^ with ek = y fc^ + M*^, 
and spin and isospin degeneration as well. One gets 

Po = ^/dA;A;2[n,^(r,/i)+n_,,(r,/i)] (45) 







1 1 



4 Idke 



7r2 







Equation (^) establishes a relation between medium density and chemical poten- 
tial and is used to determine /x for a given medium density po and temperature T. 
It is also used to test the accuracy of the numerical procedure and to determine 
the necessary size of the basis. For that aim we evaluate the baryon density for a 
homogeneous a field by means of Eq. ( ^if ) within the discrete basis and check the 
agreement with the result (^S|) obtained in the momentum basis. We increase the 
basis until sufficient agreement is reached. The result is shown in Fig. 2 (dashed 
lines). Apart from a region close to the edge of the box, which is sufficiently 
far away from the soliton, the medium density is well reproduced by the discrete 
basis with a finite number of states. 

The size of the various contributions to the solitonic baryon density ( ^2]) and 
their modification when changing the medium parameters from the vacuum to 
values close to the border of instability is illustrated in Fig. 2. The dominating 
contribution results from the valence part {a = val) of the medium contribution 
( PI ) giving rise to the bump around the center of the soliton. The residual terms 



in the medium contribution describe the polarization of the Fermi sea. Their 
contribution to the density is too small to be visible in Fig. 2. However, this con- 
tribution is located at larger distances than the valence contribution and has a 
remarkable infiuence on the soliton radius. Moreover it depends on temperature 
and density and contributes to the total baryon number. It is just this part of 
the total baryon number which is responsible for the deviation from one. The 
contribution (^3|) resulting from the polarization of the Dirac sea (dotted lines) 
modifies the density distribution but does not contribute to the baryon number 
(P^). Fig. 2 illustrates nicely the swelling of the soliton when increasing temper- 
ature and density. The mean-square radius of the soliton will systematically be 
studied in Sect. |[ 
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Fig. 2: Baryon density 
distribution of the soli- 
ton (full lines) normalized 
to baryon number B = 
An J dr r'^p{r) as a function 
of the distance r from the 
center. The broken lines 
show the reproduction of 
the medium density by the 
discrete basis. The contri- 
butions of the Dirac sea are 
given by the dotted lines. 
The right end of the curves 
indicates the size of the box 
(D = 18/M*), which is dif- 
ferent in all three cases. 



Now let us consider the consequences of introducing a chemical potential 
fJ's = + Sf^ for the soliton which is different from the /i for the homogeneous 
background field. In this case Eq. (|4^) has to be replaced by 



l/T 



p(r) 



-T 



dr (^rr 



tr[D{fi,)-'^DoilJ.y 



rr 



(46) 



In the asymptotic region far away from the center of the soliton (r ^ R) we can 
replace the quark propagator Z)(/is)~^ by the propagator Do(/^s)~^ in the homo- 
geneous field with the chemical potential for the soliton. This can be proven 
by expanding Z)(/is)~^ in Eq. (|46|) around -Do(/is)~^ (gradient expansion). As a 
result, the propagators differ only by terms which are proportional to the devia- 
tions of a and vr from their asymptotic values and by terms proportional to their 
derivatives which vanish in the asymptotic region. So we get 



l/T 

p{r:^R) = -T J dri^rr tr 



= E [M(T,/is)-M(T,/i)]$°t(r)$°(r 



rr 



(47) 



with the result that the soliton density vanishes at large distances from the center 
only if the chemical potentials /i and Hs are equal. Introducing a different chemical 
potential /ig one modifies the occupation probability for quarks in unbound states 
which contribute to observables at large distances. As a result, a finite fraction 
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of the baryon number (and of other observables as well) is uniformly spread over 
the whole space. The root mean square (r.m.s.) radius R* defined by 



R* 



is infinitely large. The occurrence of unbound quark states below critical temper- 
ature and density is a consequence of the missing confinement in the NJL model. 
The situation is different for an isolated soliton at T = 0. Here one gets the soliton 
by adding 3 quarks onto the bound valence level which does not contribute to 
the density at large distances. As soon as T > and/or go> unbound quark 
levels are involved and the lack of confinement becomes evident. 

The difference Sn between solitonic and medium chemical potential which is 
necessary to ensure B = 1 amounts to a few hundreds of keV and decreases as 
1 /D^ with increasing box radius D. The resulting solitonic density at large radii 
decreases correspondingly. It vanishes in the limit D ^oo and the effect might 
be considered as a box effect. Unfortunately that is not true. Independently 
of the box size a finite fraction of the baryon number is homogeneously spread 
outside the soliton, i. e. we have J^dr r'^p{r) ^ outside any sphere with radius 
IZ around the soliton, and the mean-squared radius (|48| ) diverges. In Ref. 0, 
the (small) deviation from the medium density outside the soliton was simply 
neglected, while it was taken into account when calculating the baryon number 
B. Similar problems will occur when calculating the moment of inertia in Sect. ^. 
That is why we tolerate a baryon number slightly different from one and do not 
introduce different chemical potentials ensuring that any local expectation value 
of the soliton vanishes asymptotically. 

There is a promising method in the literature which might be applied to fix 
the baryon number of the soliton to one without changing the chemical potential. 
In Ref. the regularized version of the baryon number in vacuum, which differs 
also from one, could be constrained after introducing the chiral radius field as an 
additional dynamical degree of freedom. In the center of the soliton, this radius 
field deviates noticeably from the chiral circle. Additionally, the constraint on 
the baryon number prevents the soliton with a space dependent radius field from 
collapsing. This method will be investigated in a forthcoming paper. 



4 Energy and radius of the soliton 

In this section we display and discuss energy, baryon number and r.m.s. radius of 
the soliton defined in Sect.^ Fig. 3 shows internal and free energy as a function 
of the medium temperature T for several densities po- While the internal energy 
represents the total energy which is necessary to generate the soliton the free 
energy disregards that part of the energy which is automatically delivered by the 
heat bath. 
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Fig. 3: Total internal {E) and free energy (F) of the soliton as a function of 
the medium temperature T for various medium densities po (pnm = 0.16 fm~^: 
normal nuclear density). The dotted lines show the contribution of the 
valence quarks to the internal energy and the dashed lines represent the energy 
3M* of 3 free constituent quarks. The calculation was performed with a 
constituent quark mass M = 420 MeV in vacuum. 

The lowest right part shows the baryon number i? as a function of T for the 
densities considered in the other parts of the figure. 



The first striking feature we want to mention is the relative independence of 
the valence quark energy on temperature and medium density, and hence on the 
constituent quark mass M*. The latter determines the depth of the well in the 
solitonic a field which binds the valence quarks. The decreasing depth at growing 
T and/or po is nearly compensated by a larger radius of the self-consistently 
determined potential well with the result that the valence level is kept at an 
almost unchanged energy of roughly 500/3 MeV . The solitonic solution of the 
equations of motion disappears if the valence level comes close to the top of 
the well in the a field. Comparing total soliton energy with the mass of 3 free 
constituent quarks we notice that the soliton energy depends more weakly on T 
and Po than the constituent quark mass. 

Comparing the free soliton energy with the results of Ref. [0 we notice differ- 
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ences up to several hundred MeV especially at larger medium density. They are to 
attribute to different assumptions concerning the occupation of the valence level 
in the homogeneous medium and to the two different chemical potentials used in 
Ref. 0. On the other hand, our baryon number which decreases with increasing 
temperature superimposes the T dependence of the soliton energy. Dividing the 
free energy by the baryon number it exhibits a slight increase with increasing 
temperature. 

The r.m.s. radii R* (48) displayed in Fig. 4 indicate a swelling of the soliton 
when temperature and density increase. At low temperature the soliton swells 
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Fig. 4: Root mean-square ra- 
dius R* of the soliton as a 
function of the medium tem- 
perature T for various values 
of the medium density po in 
units of the normal nuclear 
density pnm = 0.16fm~'^ calcu- 
lated for M = 420 MeV. 
The lower part shows the de- 
viation from Brown-Rho scal- 
ing. R* and /* are radius and 
pion decay constant at given 
medium temperature T and 
density pQ while R and fj^ de- 
note the corresponding values 
at T= po = 0- 



roughly linearly with increasing medium density. The soliton at normal nuclear 
density is by roughly 20 percent larger than in vacuum. Above 125 MeV the 
r.m.s. radius grows continuously towards the deconfinement transition. There 
are two different reasons for the modification of the soliton size in the medium: 
the increase of the radius of the self-consistent mean field and the polarization 
of the medium quarks around the soliton. The first effect is rather pronounced 
and nearly proportional to 1/M*. The polarization modifies the baryon density 
very slightly but at rather large distances from the center of the soliton and 
has therefore a noticeable influence on the mean-square radius. The effect is 
positive at lower temperatures and negative at high temperatures. It raises the 
dependence of the r.m.s. radius on the medium density and reduces its dependence 
on the temperature. A comparison with the r.m.s. radii obtained in Ref. [0 is 
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difficult because of the finite baryon density outside the sohton which inevitably 
emerges in a model with two different chemical potentials. To get a finite soliton 
radius this part was obviously ignored. In contrast to Ref. we get always a 
larger radius if the medium density increases for any temperature. 

The lower part of Fig. 4 illustrates the deviation from the Brown-Rho scaling 
which predicts R/R* ~ fl/ f-w- Apart from the immediate vicinity of the 
deconfinement transition the deviation does not exceed 10 percent. 



5 Quasi-classical energy corrections 

The soliton considered so far exhibits several undesired properties which do not al- 
low a direct comparison with the nucleon or other baryons. Due to the mean-field 
approximation the translational symmetry is violated and the soliton energy is 
contaminated by spurious center-of-mass motion. We estimate the spurious part 
of the soliton energy which is connected with quantum fluctuations around the 
artificially fixed position of the soliton by means of quasi-classical methods and 
subtract it from the total energy. The same is done for the rotational degrees of 
freedom where the restriction to hedgehog configurations introduces an alignment 
of the isospin of the soliton inducing spurious fluctuations as well. Moreover we 
introduce a collective rotation of the soliton as a whole in order to equip it with 
definite values of spin and isospin and add the corresponding rotational energy 
to the total soliton energy giving rise to a mass difference between nucleon and A 
isobar. Rotations in space and isospace are not independent of each other since 
the total isospin of the hedgehog sohton is directed opposite to its spin. Fluctua- 
tions and rotational energies in both spaces are equal and have to be considered 
only once. We perform our calculation in isospace which can simpler be treated. 

The perturbative quasi-classical approach used for the determination of spu- 
rious translational and rotational contributions to the soliton energy has been 
adopted from low-energy nuclear physics where it is denoted as pushing and 
cranking approach [0, respectively. The same correction terms can be derived 
if one includes boosted and rotating meson fields in the stationary phase approx- 
imation, which leads to the effective action of the model [||. 

First we consider fluctuations of the total soliton momentum 
P = ^ d? r ir) p q{r) which are described by the dispersion 

(((AP)^)) = ((P^)) - ((P))^ (49) 

To evaluate expectation values of P and P^ we use the regularized version of 
the extended canonical quark potential ( [29| - |38D with K,0 = v-p. It describes the 
grand canonical potential in a frame boosted with velocity v relative to the rest 
frame of the soliton. On the analogy of Eq. ( PB| ) the expectation value is given by 

dnl{T,f,;v) ----^ r^Z}(^)-i-Do(^)-i)pl =0. (50) 



dv 



-iVcTTr, 



v=0 
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It vanishes for any time-independent hamiltonian h. Squares like of a one- 
body operator (p8[) can be decomposed into a one-body operator 
[i-*^]^--^^ = / d^r g^(T') q{r) and a normal ordered two-body operator [-P^](2)- The 
expectation values of the latter can be expressed by the second derivative of the 
extended canonical potential (|29| ) and the product of two one-body expectation 
values. We get 



(51) 



v=0 



Introducing the inertial mass tensor 



M{T,fi)6ik, (52) 



v=0 

which is diagonal for spherically symmetric solitons and has identical matrix 



elements, we get for the dispersion (^9]) 



{{iAPy)) = {{[P%^))+3TM. (53) 

The minus sign in the mass definition (|5^ ) results from the anti-hermitian char- 
acter of the euclidean velocity v. Equation (^2]) defines the inertial soliton mass 
by the response of the grand canonical potential to a boost at fixed values of T 
and n- Since the variation of fl at fixed T and fi is equivalent to the variation of 
the free energy (|2^) at fixed T and baryon number B, and also equivalent to the 



variation of the internal energy (21) if B and entropy S = —dVt/dT are fixed, we 
can rewrite Eq. (|52D accordingly. However, the determination via Vl is the most 
appropriate one in our case since we have an explicit representation of the grand 
canonical potential on its variables T and /i. That is not the case for internal 



( pTD and free energy (|22|). 

In the non-relativistic limit, the dispersion (^) corresponds to the following 
energy of the translational fluctuations of the soliton 

_ mpf)) _ {{\p\))) , 3 

^trans ^j^ 2M 2 ^ ' 

While the second term describes thermal fluctuations of the soliton mass center 
in a medium with T > the first term represents the energy of the unphysical 
quantum fluctuations of the mass center which has to be eliminated from the total 
soliton energy. Fig. 5 displays this energy as a function of medium temperature 
and density. The main contribution to the center-of-mass energy stems from the 
valence quarks which are confined by the well in the mean field. The calculated 
reduction of -Ecmm with increasing temperature can be explained by the swelling 
of the soliton in accordance with Heisenberg's uncertainty principle. But there 
is only a loose relation between center-of-mass energy and r.m.s. radius (Fig. 4) 
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since the soliton radius incorporates not only the modified mean field but also 
the medium polarization. 

After an equivalent consideration for the dispersion of the isospin operator 
T = Jd^r q'^lr) tq{r), where t = T/2 denotes the single-particle isospin operator, 
we get by means of the generating function (|29|) with nO = uj-t 

^.ot - 2J " 2J + 2^ ^^^^ 

for the energy of the rotational fluctuations with the iso-rotational moment of 
inertia 



= J(T,/i)5,fc. (56) 

u}=0 



The moment of inertia is diagonal for symmetry reasons and has identical diagonal 
elements. The energy of a soliton rotating semi-classically in isospace with isospin 
quantum number T and moment of inertia is given by 

-^crank " 2 J ' ^ ^ 

The corrected energy of a soliton with isospin T and spin J = T is obtained 
by subtracting the energy of the spurious quantum fluctuations (first term in 
Eqs. (0, ^5D) and adding the cranking energy (^) to the soliton energy (^11) 

^-" = ^-^M 2J—^^J—- ^^^^ 
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The difference between the masses of A isobar (T = 3/2) and nucleon (T = l/2) 
is then given by 

AE^^ ^ E^-r - E^Z'" = ^ ■ (59) 

Evaluating the corrected sohton energy (^) the expectation value of the one- 
body operator [-P^](i) has to be calculated numerically using Eqs. (|33|^8D with 
= p^ and the regularized sea contribution (pTI). The expectation value of the 
corresponding isospin operator can analytically be determined since the single- 
particle expectation values of are the same for all quark levels independently 
of the meson fields {{a\t'^\a) = {a^\t'^\a^) = 1/2(1/2-1-1)). Hence most of the 
contributions to the expectation value cancel out each other and we get 

(([r^](i)))=iVcfi^Q + l) = ^5. (60) 

The inertial parameters Ai and JT" will be determined in the subsequent subsec- 
tions. 



5.1 Inertial soliton mass 

In this subsection, we show that the inertial mass f|52D of the soliton is identical 



with its internal energy (^) and need not be calculated separately 



M=E"' + Er'''' + 



med 



E. 



.A (61) 

Assuming spherical symmetry we get by means of the derivations in appendix B, 
which result in Eqs. (|B.11| , p.23| ), for the inertial soliton mass 



1 d'niiT,fx;v) 
3 dv-dv 



MT + M 



med 



(62) 



v=0 



with 



^ 3 dv-dv 



(63) 



v=0 



-Nr 



ds lim TTr 



1/A2 



and 

_y^mcd 



1 d'^n'^'"'^'^{T,fi;v) 
3 dv-dv 



(64) 



11=0 



-NrTTl 



\ + dl + ^7- V(a + i75T-f7r)) - AM~^d^ 



lim TTr 



^(0)"'(y + dl + ^7- V(a + i75T-f7r; 



Ao(o)-ia^ 



-AT.TTr 



A{^) ij,[{h - jj,) + -r-[h,p]) - Aoijj,) ^{ho-^) 
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with A(o)(/i) defined in Eq. ( |B.1| ). Now we exploit the invariance of the potential 
Q with respect to an arbitrary variation of the meson fields a and vr around 
the stationary point in accordance with the equation of motion (|1^). A variation 
which is in accordance with both the spherical hedgehog symmetry and the chiral 
circle respecting the boundary conditions 6a = and 6n = at small and large 
separations from the center of the soliton is given by 



6a = e r^dka 



and 



6-77 = er^dkT^ = er^dk{r 



(65) 



with an infinitesimal variation parameter e. Such a variation of the meson fields 
gives rise to the following changes 5^2™, 5^2^''''^'^ and in the mesonic and 

quark contributions to the grand canonical potential (|10]) 



m 

G 
m 



dV 



6a(r) 



m 
G 



3- / dV {a - a,) 



(66) 



AT, 



with 



5_h 

e 

6Jl_ 
e 



6(h-fif 



ds lim TTr 



-sAiO) 



1/A2 

-^TTr 



Aifi)-' 6{h - /i)2j + ^ l^hn TTr [A{Qr'6h\ 



( 6a . 6it' 
/5I — + i75r- — 



f3r-'V{a + i75T-f7r) 
ir-[h, p] = cap — i[h, r-p\ , 
, 6h^ 



2p + 17- V(cr + i75T-'r tt) — i h ,r p 



6h'^ - 2fi6h . 



(67) 
(68) 

(69) 

(70) 
(71) 



Now we introduce first 6h, 6h^ and 6{h — yu)^ and then 6VlY''^^ and ^fi'i'™"'^ into 
Eqs. ( 



6^) and get by means of the equation of motion (|Tq) and the variation 



61) of Vt"^ 



M 



ds lim TTr 



(72) 



1/A2 



-AT.TTr 



+^"0 lim^TTr[(^A(0) 



-1 
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The agreement of A4 with the internal energy (^Tj) can now be estabhshed by 
means of Eqs. ( |A.l(]| - |ArT3|) by comparing the various terms in Eq. ( [72D with the 
components (p3|-25) of the internal energy. 

The equivalence of inertial soliton mass and total mean-field energy is by far 
not trivial despite the Lorentz-invariance of the initial NJL Lagrangian. The 
approximations, the particular regularization scheme applied only on the Dirac- 
sea contribution and the presence of the medium might disturb the equivalence 
of inertial mass and total internal energy. 



5.2 Iso-rotational moment of inertia and A-nucleon mass 
splitting 

The iso-rotational moment 



1 d'nl{T,fi;Lj) 
3 duj ■ duj 



j-sea, _j_ j-med 



oj=0 



consists of the components 



sea 
A 



1 d^QY""H 



3 du; ■ dijj 



4 ^ 

^ a/3 L 



a)=0 

Rj^Sq.ep] A) 



{a\T,\P){P\n\a) 



■(a rsiP }{P \T3\a } 



with 



and 



J 



med 



A e-4M' - e-^'/^' 1 / 

■ ■ + - erfe(£„/A) - erfc(£^/A) 

y/ir £p + £q 2 ^ 

1 92fi^''^^^(T,/i;a;) 



duj-du 



oj=0 



-NrTTr 



D{f,)-%D{fir% - Do{fi)-%Do{fi)-% 



+Nc hm TTr 



4 ^ 



£a — £(3 
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(73) 



(74) 



(75) 



(76) 
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While the sea component 
bution (|76D is obtained by means of Eqs. 



has been derived in Ref. P3 



the medium contri- 
Ol , [Ol , [Algp . Since the 



single-particle hamiltonian /iq of the homogeneous medium commutes with T3 only 
diagonal elements with a° = /J*^ contribute to the corresponding terms in the iner- 
tial momenta (173, 176|). Because of lim Rj{e, e'\ A) /{e — e') = these terms vanish 

— — e'^e 

in Eq. ( [74D and the homogeneous medium does not contribute to sea component 
of the inertial moment. That is not true for a calculation in the discrete basis [EO 



with boundary conditions depending on the superspin quantum number. Here we 
have numerically to determine the inertial moment of the homogeneous medium 
and to subtract from the moment of the solitonic configuration. 

Fig. 6 illustrates the moment of inertia as a function of medium temperature 
and density. At vanishing density, the moment is nearly constant and increases 




Fig. 6: Moment of inertia 
(upper part) and A-nucleon 
mass splitting AE'aat (lower 
part) as a function of the 
medium temperature T for 
various values po of the den- 
sity and constituent quark 
mass M = 420 MeV in vac- 
uum. 

The lowest line in the lower 
part shows the splitting for 
Po = 0-5pnm obtained with the 
assumption that the valence 
level is empty in the homoge- 
neous medium. 
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remarkably only in the neighborhood of the critical temperature at 185 MeV. At 
finite medium density, the increase starts earlier. The main contribution to the 
moment of inertia comes from transition matrix elements between the valence 
level and an unoccupied level in its vicinity. At finite density, the levels around 
the valence level are sufficiently occupied by quarks representing the medium 
and the moment of inertia is remarkably reduced in comparison to the vacuum 
(Pauli blocking). The resulting moment of inertia is very small and the AN mass 
splitting (Fig. 6, lower part) is huge at low temperature and finite density. This is 
an obvious shortcoming of the model describing the medium as gas of constituent 
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quarks. In a more realistic picture, the medium quarks should be bound in 
solitons and the corresponding transition matrix elements are not blocked to 
that degree. At higher temperature, the probability of finding a hole close to 
the valence level increases. Here the blocking effect diminishes. If one keeps the 
valence level of the homogeneous medium free (n^o ^ = 0) as in Refs. ^ one gets 
big transition matrix elements to that level, and the moment of inertia is huge. 
The resulting AN splitting is negligibly small already at half of normal nuclear 
density (lowest line in Fig. 6) and further decreases if the density grows. That is 
another reason why we discarded this method of tailoring a B = l soliton. 

The quasi-classical energy corrections in Eq. (^) represent approximations to 
the first terms in an l/N^ expansion. While the quantum fluctuations behave like 
(l/Nc)^ the cranking term is proportional to l/Nc- So it is not surprising that 
the mass shift at po = obtained in our approach exhibits a similar dependence 
on T as the shift evaluated in heavy baryon chiral perturbation theory (HBxPT) 



using a 1/Nc expansion p4|. The shift is negative for nucleons and positive for A 
isobars and has the same absolute value in our approach apart from a term which 
is proportional to the deviation of the baryon number from one. The identity 
of the absolute values of the mass shifts for nucleon and A isobar is the result 
of the restriction to 2 quark flavors in contrast to the HB^PT calculation which 
includes strange quarks. At T ^ 130 MeV the AN splitting is reduced by only 
5% in comparison to 20% in Ref. p4|. Again a partial blocking of quark levels 



in the neighborhood of the valence level prevents a larger moment of inertia and 
reduces the decrease of the AN splitting at flnite temperature. 



6 Energy of the nucleon 

In Fig. 7, we display the corrected internal energy ( ^81 ) and the corresponding 
free energy in dependence on temperature and density of the medium for nucle- 
ons (T= 1/2). In the considered region, the baryon number varies between 1.2 
and 0.8 as displayed in the lower right corner of Fig. 3. To estimate the effect 
of the varying baryon number we display the energy per baryon number on the 
right panel of Fig. 7. We see that the variation in the baryon number has only 
a moderate influence on the corrected soliton energy. The behavior of the soli- 
ton energy in dependence on temperature and density differs remarkably from 
the corresponding behavior of free constituent quarks (dotted lines). While con- 
stituent quarks get lighter with increasing temperature the soliton gets heavier. 
The dependence on the medium density is weaker for solitons than for constituent 
quarks. 

The increase of the nucleon mass is mainly due to the reduction of the center- 
of-mass energy (Fig. 5) which shrinks from 350 MeV at T = to 100 MeV close to 
the critical temperature. This has to be taken into account if one compares with 
calculations which do not consider this spurious energy. A slight decrease of the 
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Fig. 7: Total corrected internal [E) and free energy (F) of the nucleon as a 
function of the medium temperature T for 3 values of the medium densities 
po and constituent quark mass M = 420 MeV in vacuum. 
The left panel shows the corrected energies for particles with the varying 
baryon numbers displayed in Fig. 3. The corresponding energies per baryon 
number are shown on the right panel. 



nucleon mass at higher temperature as e. g. observed in Ref. 25 is changed into 
an increase by means of the center-of-mass energy. Center-of-mass corrections do 
also reduce the density dependence of the nucleon mass at low temperatures. 

We should mention that the calculated nucleon mass in vacuum is by roughly 
200 MeV smaller than the experimental value. This is an obvious shortcoming of 
the simple effective model and the approximations applied in the course of the 
evaluation. For that reason the model is preferably used for the evaluation of 
the splitting between the masses of different baryons. In that sense we do not 
consider the absolute masses but their variation in dependence on temperature 
and density. Furthermore we use the experimental A-nucleon mass-splitting in 
vacuum in order to fix the only free parameter of the model - the constituent 
quark mass in vacuum - to a value of 420 MeV . 
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7 Conclusions 



We investigated the properties of a two-flavor NJL soliton which is embedded 
in a medium of constituent quarks with self-consistently determined constituent 
mass. Energy and mass of the sohton are determined in mean-field approximation 
with the restriction to hedgehog configurations and to the chiral circle. To get 
a solitonic solution of the corresponding equations of motion we have to fix the 
occupation probability of the valence level independently of the thermodynam- 
ical parameters of the medium. Otherwise the soliton dissolves below 100 MeV 
already at densities below the normal nuclear density. The expected critical val- 
ues of medium temperature and density are obtained with the assumption that 
the occupation probability of the valence level equals to one, the same value as 
assumed for the soliton in vacuum. 

Through lack of confinement the model does not allow the construction of a 
localized soliton with fixed baryon number as soon as medium temperature or 
density differ from zero. Keeping the baryonic charge confined within a finite 
radius around the soliton the baryon number of the self-consistent field config- 
uration varies between 0.8 and 1.2 in dependence on temperature and density. 
Fixing the baryon number to a definite value by means of a chemical poten- 
tial which differs from the chemical potential of the medium a part of solitonic 
baryon charge is uniformly distributed over the whole space. This is an obvious 
contradiction to the definition of a soliton. 

To remove spurious contributions to the mean-field energy and to equip the 
soliton with the quantum numbers of nucleon or A isobar we adopted the quasi- 
classical pushing and cranking approaches. The resulting energy corrections are 
determined by inertial parameters describing the response of the soliton as a 
whole with respect to a translation or rotation. We found the nontrivial result 
that the inertial mass in the medium is identical with the internal energy of the 
soliton. The rotational moment of inertia was determined numerically. 

It has turned out that the description of the medium as a non-interacting 
gas of constituent quarks moving in the solitonic mean field overestimates the 
effect of the medium on the soliton. In particular, the expected decrease of the 
AN splitting at increasing temperature and density is remarkably reduced by the 
quarks of the medium. At lower temperatures, the Pauli blocking of low lying 
quark levels by medium quarks dominates the behavior of such quantities which 
are described by transition matrix elements between different quark levels. It 
overcompensates, for instance, the influence of the swelling effect on the moment 
of inertia. Instead of increasing the moment of inertia decreases with increasing 
medium density. 

As a result of its internal structure, which is generated by a self-consistently 
determined mean field, the behavior of the soliton energy in dependence on tem- 
perature and density deviates remarkably from the corresponding behavior of the 
constituent quark mass. The scaling property between both quantities is notice- 
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ably disturbed since the influence of the changed constituent mass (depth of the 
well in the mean field) on the soliton energy is accompanied by an variation of 
the size of the well in the self-consistent mean field. 

After subtracting translational and rotational corrections the discrepancy gets 
even more pronounced since translational and rotational corrections decrease with 
increasing temperature and density. As a result the soliton mass increases with 
increasing temperature while the constituent mass decreases. 

The swelling effect of the soliton in dependence on medium temperature and 
density is well pronounced. It does not only correspond to the increase of the 
radius of the self-consistent mean field but is also related to the polarization of 
the medium in the neighborhood of the soliton. The latter intensifies the swelling 
with increasing medium density but reduces the dependence on temperature. 
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Evaluating the trace Tr of an operator 0(5,-, h) containing the differential op- 
erator dr and the time-independent operator h, which includes functional trace 
with anti-periodic boundary conditions in the euclidean time-interval [0, 1/T] and 
traces tr over Dirac and Pauli matrices, we use the representation 



Appendix 



A Operator traces 



1/T 




(A.l) 







a n=— oo 
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with the eigenvalues Sa of h and the Matsubara frequencies cu^ = (2n+l)7r T. At 
T^O, the sum over n has to be replaced by an integral 



+00 



(A.2) 



Within Schwinger's proper-time regularization scheme the regularized trace of 
the logarithm of a positively definite single-particle operator (9 at T — is given 
by 



lim TTrA \nO(dr,h) 



/ -limTTre-^^^^-'^) (A.3) 

J S T^O ^ ' 



ds 
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When calculating traces such as in Eq. ( |A.1|, |A.2|) we use the relations 



+00 



T E [ln(.;^+a^)-ln(.;^ + 6^) 



-6 + 2Tln(l+e-"/^) -2Tln(l+e- 



b/T 



2^ 



In (uj^W)-\n (oo^ + b^) 



\a\ — 161 



(A.4) 



(A.5) 



and 



+00 



1 



1 



n=— 00 
+00 



lUn + a 
du 1 



2 1 + e'^/T 
sign(a) 



271 iu + a 



(A.6) 
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Evaluating products of two thermal propagators we use 



+00 



T E T 



1 



iUn + a iujn + b 



a — b 



T 



b — a 



a — b 



+00 

E 



1 



iuJn + a iuJn + b 
1 



0(-a) -e(-6) 



sign(a) — sign(6) 
~ 2(6 -a) 



(Ai 



(A.9) 
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The following identities for traces of the operators ( |B.1| ) can be proven by means 
of the representations ( |A.1| - |A.3| ) 



ds lim TTr 



- lim TTrAlnAfO) 



Tr 
Tr 



A{fl)-'fL{h-fl) 



T d 

-^TrlnA(/x) 
2 dT 



/i 5 
29^ 



Trlny4(/i) 



lim TTr 



- lim TTr InAfu) 

2 T^O 



(A.IO) 

(A.ll) 
(A.12) 
(A.13) 



In some of the equations above we have neglected an infinitely large constant 
which vanishes if one considers the difference between two traces. 



B Evaluation of the mass tensor 



Evaluating the mass tensor (|52D we introduce the hermitian operators 

^(o)(/i) = D(o)(/i)t /^(o)(/^) = -dl + (/i(o) - /i)' (B.l) 

and 

A(o) (/i; v) = D(o) ifi; vy D^o) {f^; v) = A(o) (/x) + BIq/ - {v-p)\ (B.2) 

with D{Q){jj) from Eqs. (H, |^), D(Q){jj,; v) from Eq. ( pO]) with kC = v-p, and with 
the operators 



d_ 

dv 



v=0 



Bl 



d 



dv 



(B.3) 



(B.4) 



v=0 



which are independent of the chemical potential /i. Here we consider more general 
meson fields cr(r) and 7r(r) which are not necessarily restricted to hedgehog 
configurations and to the chiral circle. The commutator [h, p^] in Eq. ( p.3|) is 
given by the derivative of the mean field and vanishes for h = h^. Following 
Ref .|TT| we introduce the commutator representation of B^ and Bq 



B] 



(0) 



C\A(o)(0) = C\A(o)(/i) + 2/i/i(o) 



(B.5) 
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with 



a 



(B.6) 



First we treat the proper-time regularized sea contribution and notice that the 
first derivative of the exponential function is given by 



d 



sA(0;«) ^ _^ /^^g-(l-i)sA{0;-u) 



-tsA{0;v) 



(B.7) 





At = only B'^ survives in the inner bracket and can be replaced by the 
commutator ( |B.5| ). The integral is just the commutator between and 
(see e. g. appendix of Ref. p6| ) 



d 



sA{Q;v) 



11=0 



c^-sA(o) 



-tsA{Q) 



(B.8) 



The second derivative is obtained by differentiating Eq. ( |B.7|) . At v = we can 
apply Eq. (|B.8|) and get 



sA{0;v) 



v=0 



-s dt 



g-{l-t)sA(0)l ^fcg-tsA(O) 



(B.9) 



1 1 

+s /dte-(^-*)^^(°)2pye-*^^(°)-s fdte-'^'-'^'^^'^B'' [c\e-*^^(°) 



Calculating the trace of expression ( p.9|) the various terms can be rearranged and 
simplified. The integration over t becomes trivial 



Tr 



sA(0;v) 



2Tr 
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{pV + -[c\b 
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and we get 
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-Nr. / ds lim TTr 
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with the commutators 

[C\B''] 
[C\ Sol 



-sA(0) ( pipk _^ ^ 



C\B^ 



1 -sAo(O) 

2 



25^'^9^ + iY9fc(a + i75r-7r) 



:B.12) 
:B.13) 
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Notice that Tr[e~'^'^°^°-'p*p'^] vanishes because of p^p'^ = |[y4o(0), r^p'^] and the 
cychc property of the trace. 

Now we consider the medium contribution (O) to the inertial mass and find 



Tr \nA{fx;v] 



-2Tr 



v=0 



:B.14) 



To evaluate the second term we apply the commutator representation ( [B.5|) of 
the operator and get 



Tr 



A{n)-' B'A{^i)-^ B'' 



Tr 
Tr 



[C\ A{fi) + 2fih] A{fL)-^ B^ 
'A{^ir^[C\A{^,)]A{^ir'B\ 
+2/iTr[A(^)-i [C\h]A{ii)-^ B'' 
The first term in Eq. ( |B.15| ) will be treated as in |Tl| yielding 



Tr 



A{^^)-' [C\A{^^)]A{^^)-'B' 



Tr 



A{^i)-^[C\B''^ 



To reformulate the second term we rewrite the commutator 
[C\h] = - '-{r\A{^^)} + \D{^^)^r^D{^^) 
with {A, B} = AB + BA and get 



Tr 



A{i,)-'[C\h]A{^iy'B' 



-iTr 
— iTr 



+ iTr [D{^i)^) 'B''D{fi)-V 



Ai^^)-'^-{B^r^} 



Using Eqs. |B.3|) we obtain 



and 



Tr 



Tr 



Tr 



(B.15) 
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(B.18) 



(B.19) 



(B.20) 



The last term does not contribute to the trace since h commutes with ^. 
Altogether we have 

92 



Tr \nA{fi; v] 



(B.21) 
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A(/i)-i(2py + [C\B''] + 2/i [{h - ^i)S'^ + ir%,p^ 
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and 
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Tr In Ao(/i; v] 



-Tr 



Aoifi)-' ([C\ i?o1 + 2/i(/io - ) (B.22) 



v=0 



with the commutators [C*,i?^] and [C\ Bq] given in Eqs. ( [B.l"^ , |B.13| ). Finally 

we get 
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(B.23) 



v=0 
1 



-iVcTTr 



References 

[1] Nambu, Y., Jona-Lasinio, G.: Phys. Rev. 122, 345 (1961); 124, 246 (1961) 

[2] Christov, C.V., Blotz, A., Kim, H.C., Pobylitsa, P.V., Watabe, T., MeiBner, 
T., Arriola, E.R., Goeke, K.: Prog. Part. Nucl. Phys. 37, 1 (1996) 

[3] Alkofer, R.,Reinhardt, H., Weigel, H.: Phys. Rep. 265, 139 (1996) 

[4] Vogl, U. and Weise, W.: Prog. Part. Nucl. Phys. 27, 195 (1991) 

[5] Hatsuda, T., Kunihiro, T.: Phys. Rep. 247, 221 (1994) 

[6] Jaminon, M., Ripka, G., Stassart, P.: Nucl. Phys. A504, 733 (1989) 

[7] Berger, J., Christov, C.V.: Nucl. Phys. A609, 537 (1996) 

[8] Birse, M.C.: J. Phys. G 20, 1537 (1994) 

[9] Christov, C.V., Goeke, K.: Nucl. Phys. A564, 551 (1993) 

[10] Ring, P., Schuck, P.: The Nuclear Many-Body Problem. Berlin, Heidelberg, 
New York: Springer 1980 



32 



[11] Pobylitsa, P.V., Arriola, E.R., Meifiner, T., Griimmer, F., Goeke, K., Bro- 
niowski, W.: J. Phys. G18, 1455 (1992) 



[12 

[13; 

[14 
[15 

[16; 

[17; 
[is; 

[19 

[20 
[21 

[22 

[23 
[24 
[25 



Stratonovich, R.L.: Dokl. Akad. Nauk S.S.S.R 115, 1907 (1957) 

Hubbard, J.: Phys. Rev. Lett. 3, 77 (1959); Proc. Roy Soc. A276, 238 

(1963) ; Proc. Roy Soc. A277, 237 (1964); Proc. Roy Soc. A281, 401 

(1964) ; 

Sieber, P., Meissner, Th., Griimmer, F., Goeke, K.: Nucl. Phys. A547, 459 

(1992) 

Watabe, T, Toki, H.: Prog. Theor. Phys. 87, 651 (1992) 

Schwinger, J.: Phys. Rev. 82, 664 (1951) 

Eguchi, T.: Phys. Rev. D14, 2755 (1976) 

Meissner, T., Goeke, K.: Nucl. Phys. A524, 719 (1991) 

Wiinsch, R., Goeke, K., Meissner, Th.: Z. Phys. A348, 111 (1994) 

Kahana, S., Ripka, G.: Nucl. Phys. A419, 462 (1984) 

Schlicnz, H., Weigel, H., Reinhardt, H., Alkofer, R.: Phys. Lett. B315, 6 

(1993) 

Brown, G.E., Rho, M.: Phys. Rev. Lett. 66, 2720 (1991) 

Reinhardt, H.: Nucl. Phys. A503, 825 (1989) 
Bedaque, P.F.: Phys. Lett. B387, 1 (1996) 

Christov, C.V., Arriola, E.R., Goeke, K.: Nucl. Phys. A556, 641 (1993) 

Veltman, M.: Diagrammatical The Path to Feynman Diagrams, Cambridge 
Lecture Notes in Physics 4. Cambridge: Cambridge University Press 1994 



33 



